UAB-FT-708 
MPP-2012-96 



Postcards from oases in the desert: 
phenomenology of SUSY with intermediate scales 

Carla Biggio a , Lorenzo Calibbi b , Antonio Masiero c , Sudhir K. Vempati d 

a Institut de Fisica d'Altes Energies, Universitat Autbnoma do Barcelona, 08035 Bellaterra, Barcelona, Spain 

b Max-Planck-Institut fur Physik (Werner-Heisenberg-Institut), Fdhringer Ring 6, 
D-80805 Munchen, Germany 

c Dipartimento di Fisica, Universita di Padova, and INFN, Sezione di Padova, via F. Marzolo 8, 

35131 Padova, Italy 

d Centre for High Energy Physics, Indian Institute of Science, Bangalore 560 012, India 

Abstract 

The presence of new matter fields charged under the Standard Model gauge group at 
intermediate scales below the Grand Unification scale modifies the rcnormalization group 
evolution of the gauge couplings. This can in turn significantly change the running of 
the Minimal Supcrsymmetric Standard Model parameters, in particular the gaugino and 
the scalar masses. In the absence of new large Yukawa couplings we can parameterise 
all the intermediate scale models in terms of only two parameters controlling the size of 
the unified gauge coupling. As a consequence of the modified running, the low energy 
spectrum can be strongly affected with interesting phenomenological consequences. In 
particular, we show that scalar over gaugino mass ratios tend to increase and the regions 
of the parameter space with neutralino Dark Matter compatible with cosmological obser- 
vations get drastically modified. Moreover, we discuss some observables that can be used 
to test the intermediate scale physics at the LHC in a wide class of models. 

1 Introduction 

The apparent unification of the Standard Model (SM) gauge couplings can be regarded as a 
major achievement of the Minimal Supersymmetric Standard Model (MSSM). Indeed, gauge 
couplings unification represents the most convincing hint of a Grand Unified Theory (GUT) 
at very high energy scales. Typically, unification is achieved by assuming the absence of 
new physics between the electroweak (EW) (or supersymmetric (SUSY)) scale and the GUT 
scale, Mqut ~ 10 16 GeV. In fact, the presence of fields charged under the SM gauge group 
at intermediate scales below Mqut modifies the renormalization group (RG) evolution of the 



gauge couplings and in general spoils the successful gauge coupling unification. On the other 
hand, the presence of "oases" of new physics in the "big desert" between 1 TeV and Mqut 
is a natural prediction of many extensions of the MSSM. For instance, neutrino masses point 
towards a lepton number breaking scale some orders of magnitude below Mqut- The fields 
associated with such new scale can be charged under the SM gauge group, as in the case of 
the so-called type-II [lj or type-Ill [21 [3] seesaw models. Also, the dynamical generation of the 
SM flavour hierarchy typically requires heavy vectorlike quarks or Higgs fields as mediators of 
the flavour symmetry breaking [4J (for a recent discussion see [5] ) . Finally, intermediate scales 
are present in models where the breaking of the GUT symmetry to the SM one is achieved 
via intermediate steps, such as Pati-Salam models [6j or left-right symmetric models [7j. 

Gauge coupling unification can be maintained by appropriately choosing the masses of 
the new fields or by embedding them in particular sets (the simplest ones being complete 
multiplets of a GUT group, but more general choices are also possible [SJ). Nevertheless, the 
value of the unified coupling is modified by the presence of the new fields. This can in turn 
significantly change the RG running of the MSSM parameters, in particular the gaugino and 
the scalar masses (if the SUSY breaking occurs at scales higher than the intermediate scale) . 
Therefore, one can expect a potentially observable impact of the intermediate-scale physics 
on the low-energy SUSY spectrum and phenomenology. Recently, the possible consequences 
of intermediate scales have been discussed in a variety of specific models [9]-|15j. 

In this work we are going to discuss the phenomenological consequences of the intermedi- 
ate scale, due to the modified running of the MSSM parameters. We assume that the SUSY 
breaking scale is equal or larger than Mgut ; such that the running of the SUSY breaking 
masses is indeed affected by the presence of the intermediate scale. Without restricting to a 
specific model, we consider generic sets of new matter fields (i.e. chiral superfields) in vector- 
like representations of the SM gauge group, forming approximately degenerate multiplets of a 
GUT group. This allows us to highlight the common features and the possible observable con- 
sequences of this kind of models, as it ensures that gauge coupling unification is maintained 
independently of the intermediate scale and Mqut is the same as in the ordinary MSSM. 
These properties are in general not satisfied in the case of multiple step breaking of the GUT 
symmetry to the SM, i.e. in the presence of new gauge bosons (vector superfields) below 
Mqut (see however [8]). Therefore, we assume the gauge group to be 577(3) x 577(2) x U{\) 
up to the GUT scale. However, as we will see, the effects we are going to discuss mainly 
depends on the fact that, in the presence of new matter, the SM gauge couplings unify at a 
common value that is larger than in the "big desert" scenario. As a consequence, we expect 
that our findings qualitatively occur in broader classes of models, whenever the latter feature 
is realised. 

The rest of the paper is organised as follows: in section [2] we discuss the main effect, i.e. 
the modification of the running of the MSSM parameters in the presence of new fields at 
intermediate scales. In section [3] we discuss a few observables that can be used to test the 
intermediate scale at the LHC, while in sections [4] and [5] we analyse the consequences for the 
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Figure 1: Modification of the gauge couplings running in the presence of matter at the 
intermediate scale Mj = 10 13 GeV. In the left panel we take Ab = 5 (corresponding e.g. to a 
single 24 representation of SU(5)), in the right panel we take Ab = 15 (e.g. 3 x 24). 

neutralino relic density and the proton decay, respectively. Finally, in section [6] we conclude. 
Details on analytical solutions of the one loop running and a discussion on the impact of two 
loop RGEs are given in the two appendices. 



2 MSSM running with an intermediate scale 

Our starting assumption is the presence of a set of chiral superfields in complete vectorlikc 
representations of SU(5) at an intermediate scale Mj < Mgut- This choice does not spoil 
the successful, one loop, gauge coupling unification of the MSSM, as the running of the 
three couplings gets deflected in the same way. In other words, the successful prediction for 
as(Mz) (and Mgut) is not modified^ However, it is well known that the fields at Mj make 
the running above this scale "stronger" and the gauge couplings finally unify at a value ajj 
larger than in the MSSM. This effect can be seen by solving the one loop RGEs: 

1 _ 1 bf M M s 6° M GUT Ab M GUT _ 1 Ab M GUT 
au ai{M z ) 2vr n M z 2vr 11 M s 2vr n Mj ~ cfy 2vr 11 M 7 ' ^ ' 

where afj is the unified coupling in the MSSM without intermediate scale (afj ~ 1/25), 
b SM _ (41/10,-19/6,-7) and 6° = (33/5,1,-3) are respectively the SM and MSSM f3- 
function coefficients for a« (i = 1,2,3), Ab is the universal contribution of the additional 
fields at Mj and Mg is the typical low-energy SUSY scale. Ab is given by the sum of the 
Dynkin indexes of the SU(5) representations of the fields at M/j^] From Eq. ([!]), we see that, 
since Ab > for chiral superfields, the unified coupling ajj is in general larger than the MSSM 

1 This conclusion holds under the assumption that there are no large mass splittings among the fields in 
the SU(5) multiplets. 

2 For example in the SU(5) embedding of type-II seesaw [16], the new fields are in a 15 + 15 representation, 
which gives Ab = 7, while in SUSY type-Ill seesaw [171 1101 111] each copy of 24 contributes with Ab — 5. We 
remind that a copy of 5 + 5 corresponds to Ab = 1. 
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Figure 2: Contours on the plane Mj-Ab of the inverse of the unified gauge coupling 1/ajj 
(considering one loop RGEs in the left panel, two loops in the right one). The white region 
is excluded from the requirement of perturbativity of the couplings up to the GUT scale. In 
the right panel light blue corresponds to 1/ctjj > 25, blue to > 20, yellow to > 15, red to 
> 10, brown to > 5 and black to > 0. We refer to Appendix B for details on the meaning of 
A b in the two loop case. 

one, ajj > afj^ This effect is exemplified in Fig. [TJ for M/ = 10 13 GeV and Ab = 5, 15. The 
dashed lines represent the ordinary MSSM running. 

Clearly, for a given Ab, Eq. will give a lower bound on the scale Mj by requiring 
perturbativity of the gauge couplings up to the GUT scale. We can see the perturbativity 
bound in Fig. [5J where contours for 1 ja\j are plotted on the Mj-Ab planej^] The white area 
in the plot is excluded since it corresponds to 1/ajj < 0, i.e. to a Landau pole below the 
GUT scale. Looking at the left panel, where results obtained using Eq. (§ are shown, we 
see for instance that with Ab = 15 (e.g. corresponding to 3 x 24), the intermediate scale 
is constrained to be Mj > 10 12 GeV, while in case of the type-II seesaw (Ab = 7) we have 
Mi > 10 7 GeV. On the other hand, it is remarkable that M/ can be as low as the TeV scale, 
provided that Ab < 5^] We thus observe that already with this simple requirement, we can 
exclude a large part of the parameter space. 

The figure in the left panel was derived using Eq. Q but does not qualitatively change 
considering two loops RGEs, as can be seen in the right panel of Fig. [2] Notice, in particular, 
that the region far from the Landau pole is practically unaltered, while large modifications 

3 For vector superfields Ab would be negative, implying a reduction of the value of ajy and a consequent 
modification of all the effects discussed here. However, since new gauge groups are usually accompanied by new 
chiral superfields, as long as the net effect is an increment of au, the results discussed here will qualitatively 
hold. 

4 Notice however that Ab is a discrete quantity. 

5 This result is well known in the context of gauge-mediated SUSY breaking, see e.g. |18| . 
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appear for ajj > 0.2. It is evident from this plot that constraints derived at one loop are then 
conservative. Details on two loops RGEs will be given in Appendix B. 

As we will see in the following, the main effects we are going to discuss are linked to the 
larger values of ajj induced by the intermediate scale physics and therefore can be conveniently 
illustrated in terms of two additional parameters only, Mj and Ab. 

2.1 Running of gaugino masses 

Let us now move to consider the effect of intermediate-scale physics on the gaugino mass 
running. As we know, the /3-functions of the gaugino masses are related to those of the 
corresponding gauge coupling, hence the modification of the running of the gauge couplings 
above Mj will affect the running of the gaugino mass parameters, Mj (i = 1,2,3), as well. 
This effect can be easily related to the increase of the unified gauge coupling ajj . We can see 
this from the usual one loop relation among gaugino masses and gauge couplings, which is 
not modified in our scenario: 

M(/,) = M(M GUT )^, (2) 
au 

where /i is the renormalization scale. Obviously, the low-energy gaugino mass ratios (Mi : 
M2 : M3 ~ 1 : 2 : 6 in the case of gaugino mass unification) are then not modified by the 
presence of the intermediate scale. However, for given initial values M^Mqtjt)) the low- 
energy gaugino masses result smaller than in the MSSM, since ajj is larger. This effect is 
depicted in Fig. [3] for the same choices of Mj and Ab of Fig. [TJ assuming gaugino mass 
unification with Mj(Mgut) = Above the scale Mj, gaugino masses can have a very 
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Figure 3: Running of gaugino masses for Ab = 5 (left), Ab = 15 (right). Solid lines correspond 
to the running with new matter at the intermediate scale Mj = 10 13 GeV, while dashed lines 
represent the ordinary MSSM running, corresponding to the same values for Mj at low energy. 

strong running, while below the usual MSSM evolution is clearly recovered. The dashed lines 
in the plots represent the ordinary MSSM running, giving the same gaugino mass spectrum 
at low energy. It is clear from Eq. pj that, in order to obtain the same gaugino masses at low 
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energy, the starting unified mass My 2 in absence of the intermediate scale must be smaller 
than Mi/2- 

«fc-*4-*/»(l-<*£-^)- (3) 

At this stage, this might be seen as a trivial rescaling of -M]/ 2 , i.e. a larger value at Mqtjt 
than in the MSSM is required to reproduce a given gaugino spectrum in the presence of the 
intermediate scale. However, this affects non-trivially the running of other MSSM parameters, 
in particular scalar masses, as we are going to discuss in the following. 

2.2 Running of scalar masses 

In order to see the effect on the evolution of the scalar masses, let us recall the form of the 
one loop RGEs for the soft SUSY breaking mass terms. Denoting sfermion and Higgs fields 
as (f>, we schematically have: 

J t m l = ~ E C^aMMM 2 + 1^2 ( y2 ^' + A ") ■ ( 4 ) 

i 

where we can choose t = ln(/i/MGUT) and Cj(r^) is the quadratic Casimir of the representa- 
tion of the field <j). Y 2 generically denotes Yukawa couplings (with m^, we indicate that soft 
masses of different scalars can appear) and similarly A 2 schematically refers to contributions 
proportional to the A-terms. From Eq. Q, we can see the well known behaviour in the scalar 
masses evolution: in the running from Mqut to low energies the gauge part of the /3-function 
(oc ai{t)\Mi(t)\ 2 ) tends to increase the scalar mass m~L, while the terms proportional to the 
Yukawa and the trilinear couplings have the opposite effect and tend to decrease it. This 
latter effect can be however sizeable only if large third generation Yukawas and A-terms are 
involved (such as in the case of the stop and H u masses, where these terms are proportional 
to y 2 ). Therefore for what concerns 1st and 2nd generation sfermion masses, we can consider 
only the gauge term in the /3-functions and obtain simple analytical solutions of the one loop 
RGEs (see Appendix A)Q 

How is the running of scalar masses affected by the intermediate scale fields? As we have 
seen in Figs. [T] and [3j a% and Mi run to larger values above Mj. Starting with the same scalar 
and gaugino masses at Mqut and running down to the EW scale, in the presence of the 
intermediate scales, scalar masses will grow less than in the MSSM case, because of the fast 
decrease of gaugino masses shown in Fig. |3j But while considering only gaugino masses the 
MSSM spectrum could be recovered just by rescaling the GUT values M^Mqut), this is not 
any longer true if scalars are also taken into account. In other words, the intermediate scale 
has low-energy consequences in form of a distortion of the SUSY spectrum. In fact, for the 
same values of low-energy gaugino masses as in the MSSM, scalar masses feel an enhancement 

6 Later on, when we will discuss the effects of the intermediate scale on the electroweak symmetry breaking 
(EWSB) or on third generation sfermion masses, for which the Yukawa contribution cannot be neglected, we 
will solve numerically the full set of two loops RGEs, cf. Appendix B. 
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Figure 4: Ratio of the 1st or 2nd generation RH sleptons over the bino mass M\ for m^/M\ = 
0, 1, 2 at M GUT . 
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Figure 5: Ratio of the 1st or 2nd generation LH squarks over the gluino mass M3 for itIq/Ms = 
0, 1, 2 at M GUT . 

of the gauge part of the ^-function in Eq. Q in the first stage of the running, between Mqtjt 
and Mi, due to larger couplings and heavier gauginos at high energy. This means that the 
scalar masses run to larger values, or, more precisely, the presence of intermediate-scale fields 
tends to increase the scalar over gaugino mass ratios, m^/Mj. Notice that here, in order 
to restrict our analysis to model independent effects, we are assuming no large (i.e. 0(1)) 
Yukawa couplings among the MSSM fields and the new matter. In that case the scalar masses 
would receive a further negative contribution, as shown in Eq. [4j Such an assumption is also 
motivated by the protection from new large flavour violating effects: for instance in type 
II or III seesaw models the new Yukawas are typically constrained to values that have a 
negligible impact on the SUSY spectrum (< (0.1)) by the bounds on lepton flavour violating 
processes [TO], [Til ES] ■ 

The effect sketched above can be seen by looking at the RGE of such mass ratios. In 
particular, let us consider the ratio of the right-handed (RH) selectron and Bino masses, 
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Given the minus sign in the /3-function, the ratio increases in the running from high to low 
scale and, in the presence of intermediate scale physics, it increases more, due to a larger 
gauge coupling in the high-energy part of the running. This can be also shown by means of 



the one loop formulae of the Appendix A. From Eqs. (28, 30), we get: 



m%(M GVT ) 6 



+ 



56? 



bj afj , Ab a 2 (Mj) 
6? + Aba\{M s ) bl + Ab a 2 {M s ) 



(6) 



For a given low-energy value for M\(Ms) and a given high-energy starting value m ^(-/Wgut) , 
the ratio above tends to grow in the presence of the intermediate scale (i.e. increasing Ab 
and/or decreasing Mj), since afj / a 2 (Ms) strongly grows, according to Eq. ([!]). 

The effect exemplified above is general for all scalar masses. This is the main point of 
our discussion and it is represented in Figs. |4](5j where the low-energy ratios of the 1st or 
2nd generation RH sleptons, mg, over the Bino mass M\ and the 1st or 2nd generation left- 
handed (LH) squarks, rrig, over the gluino mass M3 are plotted for different GUT values of, 
respectively, m^/M\ and rUg/M^. Again, we made use of the analytical formulae given in 
the Appendix A. As in Fig. [TJ the white area corresponds to a Landau pole occurring below 
the GUT scale. The fact that low-energy mass ratios increase by introducing new matter at 
the intermediate scale (i.e. with increasing ajj) leads to potentially observable consequences 
at the LHC and affects DM phenomenology, as we will discuss in the following sections. 



2.3 Higgs soft masses and EWSB 

Apart from scalar and gaugino masses, the presence of intermediate scale clearly affects the 
running of the Higgs mass parameters fn 2 Hu and rrijj too, and thus affects the EWSB and the 
higgsino mass [i. Once correct EWSB is imposed, /j, is given (at tree level) by the well-known 
expression: 

M 2 

M 2 = -ml u + 0« u d /(tan/3) 2 ) . (7) 

The higgsino mass [i enters the neutralino mass matrix and it is thus crucial to determine 
the composition of the lightest neutralino and whether it can be a good dark matter (DM) 
candidate. 

In the usual MSSM runs to negative values at low energy due to the terms in the 
(3— function proportional to y 2 , the top Yukawa coupling. In spite of the fact that yt at 
high energy is somewhat smaller when an intermediate scale is present Q the increase of the 
scalar over gaugino mass ratio turns out to be the dominant effect. Since stop masses enter 
the terms tx y 2 in the RGEs, the net effect is that the ratio m^/Mj tends to run to even 

7 The reason is that the terms in the yt fi- function proportional to the gauge couplings increase above Mi. 
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Figure 6: Ratio ///Mi vs. Mz for a numerical scan of the parameter space with Ab=7. 



more negative values when new physics at intermediate scale is present and thus /u/Mj gets 
increased (cf. Eq. Q). This effect in shown Fig. [6] for the illustrative case of Ab=7 and 
universal boundary conditions at the GUT scale. 

The clear tendency shown in the figure generically tells us that: (i) intermediate scale 
models tend to worsen the fine-tuning problem of the MSSM, (ii) the lightest neutralino Xi 
tends to be more and more S-like. However, independently of how strong is the effect of 
the intermediate scale, this does not rule out configurations in the parameter space that give 
fi ~ Mi, i.e. with low fine-tuning and a neutralino relic density in agreement with WMAP 
observations, due to the sizable higgsino component of Xi |19j - Furthermore, non-universal 
Higgs boundary conditions at the high scale can affect this result. 



3 LHC observables 



As we discussed above, intermediate scale physics can leave a clear imprint on the low-energy 
SUSY spectrum. The question is whether this can be observed at the LHC. In order to 
isolate the effect of the intermediate scale, we need observables that are as much independent 
as possible of high-scale scalar and gaugino masses. For instance, even if the low-energy 
scalar over gaugino mass ratios grow in the presence of intermediate scale fields, as shown in 
Figs. [4] and [5j there is also a strong dependence on the high-energy initial conditions, hence 
it seems hard to disentangle the two effects. However, even at this stage, we can see that 
SUSY searches at the LHC can be affected by the intermediate scale. Let us consider, for 
instance, the first panel of Fig. [5j that corresponds to the case of vanishing squark masses at 
the GUT scale. If we do not allow m^(MGUT) < 0, such configuration clearly gives the least 
possible ratio (mg/Mj) 111111 at low energy . While in the ordinary MSSM (mg/M3) min ~ 1, 
we see that the intermediate scale can easily push the minimum ratio to values larger than 
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2. This means that the configuration M g ~ m,Q that gives the highest sensitivity in the 
LHC SUSY searches (see e.g. [20]) would not be theoretically accessible and, independently 
of the starting values for the soft masses at high energy, only the case Mg < rrig (or even 
Mg <C mg) would be possible. On the other hand, observing the case M g « m,Q at the LHC 
would give an upper bound to ajj and thus strongly constrain the Ab and Mi parameters. 

In the following we will describe other quantities that are to large extent model-independent 
and can be used to constrain the presence of new physics at intermediate scales. 



3.1 Mass Invariants 



If we assume gaugino mass unification, the gaugino and first generations sfermion masses can 
be written (at one loop) in the form: 

Mi(M s ) = A i (M s ,Ab,M I )M 1/2 (8) 
m){M s ) = m 2 f (M Gl]T )+B f -(M s ,Ab,Mj)Ml /2l (9) 



where the coefficients A{ and Bp as functions of Ab and M/, can be read in Eqs. (26) and 
(p8l [30b. 



It is clear that in the combination 

A{' s -L-^ (10) 

the explicit dependence on the GUT-scale parameters drops, if m/~, (Mqtjt) = w^(-^GUt) = 
uiq as in CMSSM-like scenarios or, more in general, in the case of GUT-symmetric initial 
conditions (a well-motivated assumption in our setup, as we are requiring unification). Notice, 
however, that Ai and Bj do not depend on Ab and Mj only, but logarithmically on the SUSY 
mass-scale Ms as well. This induces a residual dependence of the parameters Aj J on the 
initial conditions mo and Mi/ 2 that can be relevant, as we are going to show in the following. 
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In Ref. j9], it has been pointed out that mass invariants of the kind of Eq. (10) are very 
sensitive to intermediate scale fields and can be thus useful to discriminate among different 
SUSY seesaw models (see also [121 03] ) ■ Here we want to generalise that result and study 
these invariants for generic values of M/-A6. As in Ref. [9], we consider the SU(5)-inspired 
combinations: 

m?~ — m 2 ~ rn 2 ~ — m 2 ~ rn 2 ~ — m 2 - 

Al = M 2 ' Al = M 2 ' Al = M 2 ■ (U) 

Contours for these quantities on the Mj — Ab plane are shown in Fig. |7| (taking M s = 1 TeV). 

As we can see, the invariants rapidly grow for increasing ajj. This is a further consequence 

of the effect described in the previous section: m 2 ~/M 2 increases with ajj and does it more 

w 

than m^/M 2 that does not feel the contribution of SU(2) gauginos, thus grows too. 

An analogous effect occurs for the other invariants. Measurements of the SUSY spectrum at 
the LHC can be then potentially used to derive precise information on the nature of the new 
physics possibly present at intermediate scales]^] If the invariants of Eqs. ( 11 ) will be measured 



to differ significantly from the MSSM values (A^ u » 2, A% E « 21, A DL rj 17 for M s = 1 
TeV) , this might be a hint of new physics below the GUT scale (that can provide information 
on the scale and nature of it), or simply a signal that soft masses are not universal at the GUT 
scale. In principle, it should be still possible to identify the first case by reconstructing more 
than one invariant, since the three invariants exhibit robust correlations, as Fig. [7] shows. 

As in the case of atu, the above effect is strengthened at two loops, especially in the 
vicinity of the Landau pole. The typical correction is however < 10% (see Appendix B) and 
hence the plots of Fig. [7] still give a good estimate of the impact of the intermediate scale on 
these quantities. 

As we mentioned above, the numerical values of the invariants still depend on the SUSY 

£ £/ 

scale. For instance, in the range Mz < Ms < 3 TeV, we observe a variation of A^ J up 
to 60%, which of course might spoil any attempt to constrain intermediate scale with them. 
Even if it is true that once the sparticle masses will be measured, also the SUSY scale can 
be set, we look for mass invariants which are more independent on the SUSY scale. 

We consider ratios of the above invariants, i.e. simply ratios of scalar mass differences: 

2 2 

m~ - m~ 

^lu = -4 r ■ (12) 

/3/4 m 2 _ m 2_ \ ) 

f-i fa 



Compared to the quantities in Eq. (11), these ones have a much milder dependence on Ms- 
This is true in particular for A^ = (mj~ — m^)/(m^ — mi), where the variation in the 
range Mz < Ms < 3 TeV is just of a few percent. Countour plots for these invariants are 
shown in Fig. [8) Besides using the analytical one loop expressions, we numerically computed 



8 We remind that we are considering only new chiral superfields at Mi and thus we assume only the SM 
gauge group up to GUT scale. In the presence of new gauge groups, the contribution to the invariants of the 
new vector superfields can be negative, see e.g. 12, 13 , so that cancellations might occur. 
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Figure 8: Contours for the mass invariants Anf, A^r , denned in Eq. (12), for M5 = 1 TeV. 



the invariants including two loops contributions. The results are shown in Appendix B. The 
effect of the intermediate scale is again stronger than at one loop and deviates from Fig. [8] 
quite sizeably close to the Landau pole. 

We have shown that mass invariants can in principle give information and strong con- 
straints on Ab and Mj. They might even exclude, if observed to be close to the CMSSM 
predictions, the presence of fields charged under the SM gauge group at scales below Mqtjt- 
However, it is difficult to say whether all the SUSY masses entering these invariants can be 
measured with sufficient precision at the LHC. Clearly this question depends on the actual 
mass-scale of the SUSY particles, as well as on some features of the spectrum, and is beyond 
the purposes of the present discussion. Let us only notice that, if the slepton masses can be 
reconstructed from cascade decays of heavier particles (see the next section), the invariants 
involving squark and slepton mass differences should not be difficult to obtain, given the 
typical hierarchy between coloured and uncoloured sfermions. Moreover, Figs. [7] and [8] show 
that, even in presence of uncertainties as large as 10%, the invariants can provide very useful 
information on the intermediate scale and discriminate among different scenarios. It might 
be much more difficult to resolve experimentally squark mass differences like mj~ — mj~ . How- 
ever, we notice that intermediate scale physics might help in this sense. In fact, the relative 
mass-splitting (??Iq — m^j/mg (that can at most be ~ 5% in the CMSSM) increases with 
au too and can reach values larger than 10%. 



3.2 Kinematic Edges in Cascade Decays 

Let us now investigate whether the distortion of the spectrum can be studied by means 
of kinematic observables potentially measurable at the LHC experiments. We consider the 
typical cascade decay depicted in Fig. [9j The intermediate particles are real if the following 
condition is fulfilled: 

rev, > m^o > m~ > m. o . (13) 
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Figure 9: An example of a cascade decay. 
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Figure 10: Value of the kinematic edge rnj^ in the cascade decay of Fig. |9| on the plane 
{me;, Mi) for Ab = 0,7, 15 at Mj = 10 13 GeV, respectively. 



Consequently, the invariant-mass distributions of the outgoing SM particles (jets and isolated 
leptons) exhibit sharp kinematic end-points [21]. Notice that, depending on the spectrum, 
zero, one or two sharp edges can be present. Indeed, if for example both and satisfy 
the above inequality, then two edges could be observed, while only one will be there if only 
one of the two -typically m, - satisfies it. The position of the end-points of the distributions 
can be expressed as a function of the SUSY masses: 
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(14) 
(15) 
(16) 



and can therefore be used to reconstruct the SUSY spectrum [21]. As any combination of 
SUSY masses, these observables are also modified in the presence of the intermediate scale. 

As an example, in Fig. 



10 



we plot contours for mf^ x in the plane of the physical masses 
(m^,, Mi « m x°) f° r three different choices of Ab and a fixed Mj. As expected, the value 
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Figure 11: Maximum number of edges in the cascade decay of Fig. [9] on the plane (Mj, Ab). 



of the edge changes with the intermediate scale. Again, if it is possible to independently 
measure the RH slepton and the neutralino masses (e.g. from m^ ax and m^ ax together with 
jets and missing Et distributions) and the edge of meej, this could provide a further way to 
test intermediate scale physics. 

As we have mentioned above, depending on the spectrum (in particular whether the 
sleptons are lighter than the second neutralino) the number of edges can vary. Indeed, as we 
have extensively discussed, the effect of the intermediate scale is precisely that of increasing 



the ratio of the scalar over gaugino masses and thus making the condition of Eq. ( 13 ) more 



difficult to be satisfied. In order to illustrate this point, we have obtained the maximum 
number of possible edges in the mu and m« invariant mass distributions as a function of 
(Mi, Ab). This quantity is independent of the details of the spectrum. The result is shown in 
Fig. [TT] As we can see, the effect of the intermediate scale can be such that the cascade decay 



of Fig. p| is never kinematically allowed (in other words > rn^o always) or just for the 

lighter sleptons for any choice of sfermion and gaugino masses at high energy. It is then clear 
that, if for instance two clear edges will be observed in di-electron or di-muon distributions, 
we will be able to exclude a large portion of the (Mj,Ab) parameter space. On the other 
hand, if no edges are observed at all, intermediate scale physics would remain unconstrained 
as it is always possible to choose high-energy initial conditions such that m 



X-2 



4 Neutralino Dark Matter 



The modification of the SUSY spectrum described in the previous sections can destabilise the 
regions of the parameter space where precise relations among the parameters are required in 
order to fulfill the WMAP bound on DM relic abundance. 

It is well known that the lightest neutralino of the MSSM, Xi, typically B-like (i.e. m^o ~ 
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Figure 12: Left: region on the (Mj, Ab) plane where the correct relic density for x± is obtained 
via f coannihilation (blue points corresponds to tan/3 = 10, red to tan/3 = 30). Right: the 
same for the A-funnel region. The grey dots in the background mark the region consistent 
with the perturbativity bounds. 



Ml), is overproduced in the early universe, unless the neutralino (co) annihilation cross-section 
is enhanced by particular conditions. Such conditions define few regions of the parameter 
space where the WMAP bound is satisfied: (i) the f coannihilation region, where the correct 
relic density is achieved thanks to an efficient f~xl coannihilation, which requires ~ m^o 
|22j : (ii) the "focus-point" region, where the Higgsino-component of Xi is sizable, i.e. fj, ~ 
Mi [19]; (iii) the A-funnel region, where the neutralino annihilation is enhanced by a resonant 
s-channel CP-odd Higgs exchange, if itia — 2 x m^o [23] . 

Let us consider first the f coannihilation region. In the CMSSM, such a region is usually 
a thin strip which runs along the border of a wide region of the parameter space excluded 
because it gives a f as the lightest SUSY particle (LSP) (m^ < m^o). As we have seen, the 
modification of the running due to the additional fields tends to increase the scalar masses 
compared to the gauginos. The lightest f mass is approximately given (under the conditions 
m~ 3> m| , fj, m T tan /3) by: 

2 _ 2 (nm T tan f3) 2 

m h ~ m f R 2 ' ( 17 

As we have seen, the ratio m^ R /Ml can be strongly increased by the intermediate scale. In 
particular, it can become larger than one, even for vanishing rnf R at Mqtjt (see Fig. 
As a consequence, the region with a f LSP tends to be reduced and can even disappear. 



In fact, even though the second term of Eq. (17) tends to decrease m| , the intermediate 



scale can easily make the condition < m^o impossible to obtain for any choice of the 
SUSY parameters at the GUT scale. This has been observed in Ref. |24j for a qualitatively 
similar scenario, where the SU(5) RG evolution from a universality scale above Mqtjt have 
been considered. For the reasons explained above, this is a general consequence of models 
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with intermediate scales]^] Clearly the coannihilation strip gets modified as well: it can be 
distorted or it can even disappear (in the case of a too large increase of mf R /M\ such that 
the condition ~ m x\ cann °t be achieved anymore). This has been discussed in Ref. |25j 
again in the context of SU(5) RG running. In Ref. |10j the same effect has been studied in a 
well-motivated case of intermediate-scale physics: a type I+III seesaw model, achieved with 
a single SU(5) adjoint representation, 24 (corresponding to Ab = 5). It was shown that, for 
certain choices of the parameters, the resulting coannihilation region is bounded from above 
(i.e. neutralino DM is only possible in a limited range of the neutralino mass). 

Now let us generalise these observations to generic sets of intermediate-scale fields. In the 



left panel of Fig. 12, we show the regions on the plane (Mi, Ab) where the coannihilation can 
take place. The plot was done solving the RGEs numerically at two loops with CMSSM-like 
boundary conditions. The blue points correspond to tan /3 = 10, the red ones to tan (3 = 30 
(the grey dots in the background mark the region consistent with the perturbativity bounds) . 
Comparing Fig. 12 with Fig. [2| we see that already values of ajj as large as 1/15 -f- 1/10 are 
enough to off-set completely the f coannihilation condition, since results consistently 
larger than m^o everywhere in the parameter space. The effect can be partially relaxed 
in the case of very large tan f3 (and large A-terms) but intermediate scales configurations 
corresponding to au ^ O(0. 1) will still make the coannihilation region disappear. Dropping 
the universality assumption at Mgtjt> one can still find corners of the parameter space where 
the coannihilation is possible (e.g. with m# d (MGUT) 3> mf R and large tan/3). 

The A-funnel region can face a similar fate. In fact, the condition m>A ~ 2 x m^o can be 
made theoretically unaccessible. At tree level the CP-odd Higgs mass is approximately 

rnA-mjj -m 2 Hu , (18) 



As we have seen in section 2.3 
increased too. The result is depicted in the right panel of Fig. 12 , where we show the region 



the ratios \m H \/M\ grow with au, so that tua/Mi gets 



of the plane (Mi, Ab) where the A-funnel can be obtained for tan (3 = 45. Again, if au is too 
large the condition rriA — 2 x m^o can be never realised. We notice, however, that the funnel 

-VI 

region can be restored by choosing proper non-universal values for mn u . at the GUT scale, 
as rriA and [i then become free parameter. 



Finally, let us comment about the focus point region. As we discussed in section 2.3 
the ratio fi/Mi tends to increase as well, however Fig. [6] shows that it is always possible 
to find configurations with [i ~ M\, such that the Higgsino component of Xi is sufficiently 
large to give a sizeable annihilation cross-section. The focus point region is therefore the only 
DM branch which is not destabilised by the intermediate scale, if CMSSM-like boundary 
conditions are assumed. Let us remark that this true under our assumption that the new 
fields do not have large Yukawa couplings with the MSSM fields. On the contrary, if this 
occurs, the focus point region is drastically affected and can even disappear [25J. 



J See for instance Refs. [TTI [T3l PI IT5] . 
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5 Proton decay 



In SUSY GUTs proton decay is typically induced by dimension five operators generated by 
the exchange of coloured Higgs triplets. Once a mechanism to suppress it is added, the 
model is safe, since the contribution of the dimension six operators from the gauge bosons 
associated to the unified gauge group is typically below the current experimental bound set 
by SuperKamiokande: r(p — > e + 7r°) > 1.29 x 10 34 yrs at 90% confidence level [26] . However, 
when intermediate scale physics is present, the enhancement of the unified gauge coupling 
increases the proton decay rate induced by the GUT gauge bosons to values close to the 
current bounds [151 [27] . This can thus be used to set further constraints on the intermediate 
scale. 

The partial decay width of the dominant decay mode is given by [2 
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[{A$f + {Af)\l + \V ud \ 2 ?]. (19) 



Here an = —0.0112, D = 0.80 and F = 0.47 are parameters related to the hadronic matrix 
elements and Mx is the mass of the gauge bosons of the unified gauge group, which we 

(i) 

will take equal to Mqijt in our evaluations. A R are the renormalization factors given by 
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(20) 



with 73 = 3, 72 
the proton life time: 



3 7 (D 
2 ' a 



11 7 f) 
30' '1 



U . The above expressions can be used to estimate 



r(p — > 7r°e + ) 



2.2 x 10 36 yrs 



1/25 



M x 



From this expression we see that the experimental bound r(p 
satisfied unless the unified coupling is very large ajj w 0.5. 



(21) 

2 x 10 16 GeV/ ' 

e+vr ) > 1.29 x 10 34 yrs is 



In Fig. 13 we use Eqs. (19, 20) to plot the proton lifetime on the Mj-Ab plane. The grey 
area corresponds to the current exclusion limit. As expected from the above estimate, we see 
that the present experimental bound only excludes a small portion of the parameter space, 
close to the Landau pole. On the other hand, 0(1) variations of the GUT vector bosons 
mass Mx drastically affect the predicted r(p — > vr°e + ), as the estimate of Eq. (21) shows. 
Therefore, a proton life time in the reach of the future experimental sensitivity can be easily 
obtained for large ranges of the parameters if new fields at intermediate scales are present. 
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Figure 13: Proton lifetime with gauge boson mass Mx = Mqut — 2 • 10 16 GeV. The grey 
region corresponds to the current exclusion limit. 



6 Conclusions 

In this paper we have studied the phenomenological consequences of the presence of new 
physics at a scale intermediate between the EW and the GUT scale in a SUSY theory. 
We have assumed that the new physics consists only of chiral superfields in complete GUT 
multiplets, such to maintain gauge coupling unification. 

The main effect, which drives all the others, is the increment of the value of the unified 
gauge coupling. The simple requirement that it remains perturbative up to the GUT scale is 
already enough to exclude a large portion of the parameter space, as we have shown in Fig. [2} 
As a consequence of this increase, the entire low energy spectrum is modified with respect 
to the MSSM one and, in particular, the ratio of scalar over gaugino masses is enhanced. 
This has interesting consequences both for what concerns the collider phenomenology and 
the neutralino dark matter. 

We have analysed two main sets of collider observables that can give hints of the presence 
of the intermediate scale or, on the contrary, can be used to constrain it: the mass invariants 



defined in Eqs. (10)-(12) and the edges in cascade decays. By measuring the sparticle masses 
at the LHC and building the invariants, one can in principle disentangle if we are in the 
presence of a CMSSM-like spectrum or if intermediate scales are present or if high energy 
boundary conditions are not universal. The same can be done if independent measurements 
of sparticle masses and the position of the edges in the invariant masses in cascade decays 
are available. Still, if this is not the case, we have shown that the maximum number of 
edges in these decays can give uncontroversial information that can be used to constrain the 
intermediate scale physics. 
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On the other hand, the generic increase of the ratio of scalar over gaugino masses tends 
to destabilise the regions of the parameters space where the correct DM relic density is 
obtained thanks to an efficient (co)annihilation of the lightest neutralino. This is the case 
for the f-coannihilation region or the A-funnel one: we have shown that the presence of the 
intermediate scale can render impossible to realise the precise relations among the masses 
of the involved particles necessary to enhance the (co) annihilation cross section. On the 
contrary, in spite of the presence of the intermediate scale, we found that it is always possible 
to find regions in the parameter space where the Higgsino-component of the neutralino is 
enough to increase the annihilation cross section and obtain the correct relic density ("focus- 
point" region). 

Finally, we have observed that the increment of the unified gauge coupling can reduce 
the proton lifetime if the decay is driven by the GUT gauge bosons. We have shown that, 
for gauge boson masses equal to the GUT scale, the actual bound can be used to exclude a 
small part of the parameter space. 
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A Analytical one-loop formulae 

In this section we collect analytical formulae that can be used to illustrate the effect of 
the intermediate scale on the renormalization group evolution of the MSSM parameters (see 
Sections [2]) . Aiming at compact and simply readable expressions, we consider, besides gauge 
couplings and gaugino masses, only first generations sfermion masses, for which Yukawa 
couplings can be neglected. Even though it is possible to obtain analytical solutions of the 
stop and Higgs masses as well, they result quite involved, hence we prefer to study these 
parameters numerically. 

We first write the solution of the (one loop) gauge coupling RGEs in presence of fields 
at Mj giving a contribution Ab = Y2i n i (with rii being the Dynkin index of the SU(5) 
representation of the i-th. field) to the /3-function coefficients: 




fi > Mi 



(22) 



where b\ are the ordinary MSSM coefficients 




(23) 
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and 
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Analogously for the gaugino masses we have: 
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In particular, we have: 
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For the scalar masses, if we neglect Yukawa and A-term contributions, the solution in 
general looks like: 



m?Jfj,) = m 2 JM GVT ) + Ki{n) + K 2 {n) + K 3 (fi) 



with 
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(28) 



(29) 



and for the quadratic Casimirs C{ we have: C\ = 3/5 y?, with y^ the hypercharge of 0, 
C 2 = 3/4 (0) for SU(2) doublets (singlets), C 3 = 4/3 (0) for SU(3) triplets (singlets). 



Using Eqs. (22, 26), we can easily solve the integrals and we get for /i < Mj (in particular 
fi = Ms): 
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Figure 14: Effect of two loop running on 1/ajj. The red points correspond to Ab = 5, the 
green ones to Ab = 7 and the blue ones to Ab = 15. In the three cases, the crosses correspond 
to one loop results, the diagonal crosses correspond to the two loop result with 77.24 = 1 (red), 
?ii5 = 1 (green) and 7124 = 3 (blue) and the dots show the two loop result considering only 
an appropriate number of 5 + 5. 



B Effect of two-loop RGEs 



At two loops, the RGE for the gauge couplings read: 

% = Tsfc b > + tL fe # - E crrf) , (3D 



dr l (4tt) 2 ' (4tt) 

where the two loop /3-function coefficients for the MSSM can be found, for instance, in 
Ref. [28]. Unlike in the one loop case, the new contribution depends in principle on the exact 
field content at the intermediate scale. For instance, we have: 

A6 (2) = (32) 

/ I5 n s + f n io + + f "24 §n 5 + fnio + ^nis + 15n 24 ffn 5 + f n w + ^§n 15 + f n 24 \ 

|n 5 + gnio + ^^15 + 5n 2 4 7n 5 + 21nio + 69rii 5 + 45n 2 4 16nio + 16ni 5 + 16n 2 4 

\ K n 5 + l n io + ff«i5 + f «24 6ni + 6ni 5 + 6n 24 f n 5 + 34m + ^fVs + irn 24 / 

where (x = 5, 10, 15, 24) denotes the number of 5 + 5, 10 + 10, 15 + 15 and 24 represen- 
tations at Mi, respectively. 



The effect of two loop RGEs is illustrated in Fig. 14, where the cases Ab = 5 (red), Ab = 7 
(green) and Ab = 15 (blue) are compared. In the three cases, the crosses correspond to \jau 
vs. Mi at one loop. The diagonal crosses correspond to the two loop result with the following 
field contents: 7124 = 1 (red), 77-15 = 1 (green) and 77,24 = 3 (blue). Finally, the dots show the 
two loop result considering only an appropriate number of 5 + 5, i.e. 715 = Ab. As we can 
see, the difference between "one loop equivalent" field contents (e.g. 1 x 24 and 5 x (5 + 5)) 
is appreciable only close to the Landau pole. Therefore, the two loop effects are to a very 
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Figure 15: Effect of two loop running on the invariants A® u , A® E , A DL for Ab = 5 and 
Ms = 1 TeV. The red points correspond to one loop RGEs, while the blue ones to two loops. 
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Figure 16: Effect of two loop running on the invariants Anf and A^ for Ab = 5 with 
variable Mg. The red points correspond to one loop RGEs, while the blue ones to two loops. 



good approximation independent of the exact field content and we are going to implement 
the two loop RGEs taking always = Ab. 

In the second panel of Fig. [2j 1/ctu is shown in the Ab — Mj plane. We can see that atjj 
increases with respect to the one loop case. Therefore, we expect all the effects discussed 



in the paper to be further enhanced at two loops. For illustration, we show in Fig. 15 



_ a 
7/' 



comparison between the one loop and the two loops results for the mass invariants A 
discussed in section [3] for Ab = 5 and Ms = 1 TeV. The two loop correction is mild (typically 
< 10%) and grows close the Landau pole as expected. The two loop effect can be much 
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more dramatic in the case of the quantities A^-j? defined in Eq. (12). This is depicted in 



Fig. 16 again for Ab = 5 and a variation of the SUSY scale Ms < 2 TeV. While Ms has 
no strong impact on the mass invariants (contrary to A- J ), the two loop contribution can 
largely increase the effect of the intermediate scale on the quantities A^^. Therefore the one 
loop expressions used for Fig. [8] just provide a (conservative) estimate in the regimes where 
the effect is mild and two loop RGEs should be taken into account for precise quantitative 
studies of this kind of observables. 
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